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Abstract — The mean free path (MFP) in silicon nanowire field 
effect transistors limited by surface roughness scattering (SRS) is 
calculated with the non-perturbative approach utilizing the non- 
equilibrium Green's function method. The entrance scattering 
effect associated with finiteness of the channel length is identified 
and a method to eliminate it in the calculation of the MFP is 
developed. The behavior of the MFP with respect to channel 
length (L), channel width (W), and the root-mean-square (RMS) 
of the surface roughness is investigated extensively. Our major 
findings are that the single parameter, RMS/W, can be used as 
a good measure for the strength of the SRS effects and that 
the overall characteristics of the MFP are determined by the 
parameter. In particular, the MFP exponentially decreases with 
the increase of RMS/W and the MFP versus the gate electric field 
shows a distinctively different behavior depending on whether the 
strength of the SRS effects measured by RMS/W is smaller or 
greater than 0.06. 

Index Terms — nanowire, MOSFET, nanowire field effect tran- 
sistor, quantum transport, mobility, mean free path, surface 
roughness, non-equilibrium Green's function. 

I. INTRODUCTION 

The Si nanowire field effect transistors (SNWFETs) with 
all around gates have recently emerged as a promising device 
to replace the conventional planar metal oxide semiconductor 
field effect transistors (MOSFETs) as the latter face the physi- 
cal and technological challenges in their scaling-down 0~|-|[3]]. 
The SNWFETs have the merit of good electrostatic control 
and hence can deliver high on-state currents. As their overall 
performance is determined by the aspect ratio of the channel 
length (L) and channel width (W) H, it may be desirable 
to utilize thinner and thinner nanowires. However, rough 
interfaces between Si and dielectric are unavoidably generated 
during the fabrication processes, and currents become limited 
by surface roughness scattering (SRS) as the nanowire cross- 
sectional size is reduced below 5 nm |5), 0. 

The SRS effects in SNWFETs have been a subject of 
intense theoretical and computational studies. Most of the 
studies have been performed with the semiclassical approaches 
such as the Kubo-Greenwood approach Q-|9), the multi- 
subband Monte Carlo technique ifTOl . ifTTl . and the direct 
solution of the one-dimensional Boltzmann transport equa- 
tion IT2l . |fl3l . The Schrodinger-Poisson equations are self- 
consistently solved in the approaches assuming an infinitely 
extended homogeneous nanowire, before the nanowire sub- 



Fig. 1: (a) Schematic diagram of the silicon nanowire field effect transistor 
with all-around gates, (b) Cross section view at the center of the channel with 
rough interfaces. 



band information is extracted and fed to the semiclassical 
procedures. There have also been a few studies based the non- 
equilibrium Green's function (NEGF) method OH, lfl4l-|fT9l, 
where coherent electron transport across the channel region 
with rough Si/dielectric interfaces is considered fully quantum- 
mechanically. In the non-perturbative approach, the screening 
effect, electron phase coherence, and the mode mixing effect, 
among others, are naturally considered whereas these are hard 
to be implemented directly or disregarded entirely in the 
semiclassical approaches. 

Most of the NEGF studies have so far focused on the low 
field electron mobility in nanowires. However, the mobility in 
the device with a finite channel length is not a proper measure 
for the device performance, because the mobility increases as 
the channel length increases whereas there should be more 
detrimental scattering. As Gnani has pointed out, the mean 
free path (MFP) should be a better indicator of the device 
performance ||201 . 

In this work, using the non-perturbative NEGF approach, 
we have for the first time investigated the SR-limited MFP 
in SNWFETs. Novelty and contributions of this work are: 1) 
we have devised a systematic way to calculate the MFP in 
the device with a finite channel length, 2) we have identified 
the entrance scattering effect which is unavoidable in the 
device with a source/channel/drain structure, estimated its 



strength, and eliminated it from the MFP calculation, and 3) 
we have thoroughly investigated the behavior of the MFP 
with a particular emphasis on its dependence on the strength 
of the SRS effects. These will be discussed in detail in the 
following sections. 

II. Methodology 

A. Simulated devices and Hamiltonian 

The simulated devices in this work are three-dimensional 
(3D) rectangular silicon nanowire FETs with all-around gates 
as shown in Fig. Q] The source and drain regions are heavily 
n-doped with doping concentration of 10 20 cm -3 whereas the 
channel is intrinsic. The SiC>2 oxide thickness is assumed to 
be 1 nm. To describe the electron transport in the conduction 
band, the effective mass Hamiltonian was used, with the 
effective masses being scaled according to the tight-binding 
calculation (2TJ. All the six equivalent valleys are included in 
the simulations. 

B. Surface roughness generation 

The surface roughness (SR) was realized at four Si/Si02 
interfaces in the channel region only. SR was generated ac- 
cording to the exponentially decaying autocorrelation function 
expressed as 



C(f) = A^exp(-V2r/i r 



(1) 



where A m is the root-mean square (RMS) fluctuation of the 
roughness, L m is the correlation length, and r = \r\ is the two- 
dimensional distance between two points on the interface. By 
Fourier- transforming the autocorrelation function to obtain the 
power density spectrum (PDS), multiplying a random phase 
factor to the square root value of PDS, and inverse Fourier- 
transforming the randomized PDS, instances of rough S1/S1O2 
interfaces were randomly generated lfl"&l . El . 

C. NEGF method 

We have employed the NEGF method to calculate the charge 
density and the current in the devices with rough surfaces in 
the channel region. Self-consistent calculations coupling the 
charge density from the NEGF equations and the potential 
from Poisson's equation were carried out. The methodology 
is well established in the literature lTT2l . J23l . 

The NEGF procedure with the rough surfaces is the same 
as the standard non-diffusive calculation procedure. There is 
no need to calculate the less-than or greater-than Green's 
functions, because the surface roughness scattering (SRS) 
is not treated by including the self-energy in the Green's 
function. Rather, the conduction band edge Ec{r) at the 
Si/Si02 interfaces is spatially varied in accordance with the 
generated rough surfaces. In other words, in the Hamiltonian 



H = H Q + E c (r) - q Q(l 



(2) 
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Fig. 2: (a) The effective transmission coefficient Tes f° r different wire 
widths, (b) The conduction-band profiles (1st subband) for different gate 
voltages and their corresponding transmission functions. Device dimensions 
are W = 5 nm and L = 40 nm. 



Ec{r) assumes the value of of 0.56 eV or 3.63 eV, if f 
corresponds to silicon or oxide points, respectively. In the 
Poisson's equation, 



^ • (e(r)^>(f)) = gt,(n(f) - N D ), 



(3) 



where Hq is the kinetic part of the Hamiltonian with the 
inverse effective mass tensor which is not altered by SR, 



e(f) is also spatially varying in the same manner as Ec(r)- 

In this work, the coupled mode space approach was used 
EH, ll25l with the number of modes ranging from 30 to 50 
depending on the cross-sectional width of nanowire. Including 
higher modes little changed the outcome of the calculations, 
which validates that the used number of modes is sufficient. 



III. ENTRANCE SCATTERING EFFECT 

In the nanowire FET in Fig.Q] electrons experience quantum 
mechanical scattering by the channel potential barrier as they 
enter the channel region. The channel potential barrier is 
present due to the fact that the source/drain and channel 
regions are doped with different doping concentrations. The 
entrance scattering effect (ESE) takes place even in the ballis- 
ticity condition (no rough surfaces in the channel region) and 
should be separated out from the calculated mean free path. 
The method to isolate the pure contributions from the channel 
region is elaborated in section IV. 

In the nanowire device under consideration, the ballistic 
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Fig. 3: The drain current versus channel length L for different gate voltages. 
The currents are normalized by the current of L = 55 nm. Dashed lines are 
the fit to the function of the form A/(B + L). The curves for Vq — Vt = 
-0.25 V and 0.15 V are shifted by -0.2 and 0.1, respectively, for a clearer 
view. Device width and SR parameters are W = 3 nm, A m = 0.3 nm, and 
L m = 1.0 nm. 



current If, can be expressed as, at low V 1 



DS, 



2ggVbs 



E^-i( 



Ep — Ey 
k B T 



), 



(4) 



where T m {x) is the m-th order modified Fermi integral, Ek 
is the energy level of subband k in the channel region, Ep 
is the Fermi energy, and 71- is the transmission probability in 
subband k that accounts for backscattering due to ESE. If the 
source/drain regions are removed and the same cross-section 
of the channel region is repeated indefinitely, there is no ESE 
and 71 = 1 in (0). Let us denote the ballistic current in the 
infinitely long homogeneous nanowire as I^. We write 



lb = TEshoo , 



(5) 



where the effective transmission coefficient Tes is introduced 
as a measure for the degree of ESE in the nanowire device. 

Fig. |2] (a) shows Tes versus Vq for various W's, where it 
is shown that Tes is as low as about 0.8 and its dependence 
on Vg is not sensitive to W . As shown in Fig. [2] (b), the 
channel potential barrier is so high in the subthreshold 
regime that an electron is either completely reflected or 
completely transmitted depending on whether its energy is 
below or above the potential barrier, except for a very narrow 
window of energy located just above the potential barrier 
where it experiences the quantum-mechanical backscattering. 
So Tes is close to 1 in the subthreshold regime. Around 
Vt, the potential barrier is lowered and at the same time 
its shape near the ends of the channel becomes sharper 
so that the quantum mechanical backscattering starts to 
occur seriously. In other words, the above-mentioned energy 
window broadens, resulting in a sharp decrease of Tes- If 
the gate voltage is further increased, Tes reaches a minima 
and then slightly increases, because the channel potential 
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Fig. 4: The mean free path \,sr for different channel lengths. Asr as 
defined in 41 It is also displayed. Device width and SR parameters are W — 
5 nm, A m = 0.2 nm, and L m = 1.0 nm. 



barrier is sufficiently lowered to result in lesser backscattering. 

IV. MEAN FREE PATH 

In devices with finite channel length, the effective mobility 
is no longer a good indicator of device performance, as it tends 
to zero as the ballistical condition is approached whereas the 
ON current increases lEUI . Instead, the mean free path (MFP) 
that is directly related to electron transmission in the device 
is better suited as the indicator of device performance. In this 
work, the SR-limited MFP is calculated as follows. 

At a low drain bias, we relate I with /boo via 



I = TIi 



boo i 



(6) 



where T is the overall transmission coefficient in the device. 
T consists of contributions from the entrance scattering (Tes) 
and the surface roughness scattering (Tsr)- In an attempt to 
isolate Tsr, we write 

1-T I -Tes , I -Tsr 



T 



Tes 



Tsr 



(7) 



which is written based on the simple picture that the two 
sources of scattering are connected serially ll26l . Using the 
defining relationship for the SR-limited MFP Asr, 

Asr 



we finally obtain 



Tsr 



Asr. — 



Asr + L 



L 
T-i - T -1 ' 

' ' ES 



(8) 



(9) 



Notice that in defining T in ©, hoa is used instead of lb- If 
the relationship I = T'h is used instead of (O, it implies that 
ESE is disregarded, and we may write 
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Fig. 5: The mean free path Xsr for different wire widths. The channel 
lengths were fixed at L = 100 nm for all W's. SR parameters are A m = 0.3 
nm and L m = 1.0 nm. Mean value and sample standard deviation (dispersion 
bars) are reported. 



The MFP Ksr of the above equation is related to Xsr as: 

1 
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SR- 
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That is, the MFP is multiplied by the factor 1/Tes, which 
leads to an unphysical peak as demonstrated below. 

Before proceeding further, we first show /as a function of 
L in Fig. [3] which verifies that diffusive transport expressed 
by Eqs. (__) or (fTOb actually takes place due to SR. Note that 
we have deliberately chosen the device width and the SR 
parameters in the figure so that resultant MFP is as small as 
about 10 nm (consult Fig. EJa)). 

Fig. __] shows thus-calculated \ SR for SNWFETs of W = 5 
nm. For L varying from 60 nm to 250 nm, Xsr seems to 
depend on L weakly. The MFP calculated by using ( fTTb is 
also displayed in the figure: Asr shows a pronounced peak 
around Vr when there is no physical reason for the feature. 
This supports that Xsr is a more natural choice for the MFP 
than Asr- In the Xsr versus Vg curve of Fig. [4] however, 
there is still some hint of a peak around Vr- We believe that 
this is an artifact which is originated from the imperfect nature 
in the process of isolating the pure SR contribution in the 
channel region ((_7]l). As a matter of fact, ([_} should hold better 
for greater L because the assumption that the two sources of 
scattering (ES and SR) are incoherently connected becomes 
more valid. 

Fig. [f_ shows the MFP for different W's when L is fixed 
at 100 nm. The MFP decreases as W becomes smaller, as 
expected, but the above-mentioned peaky feature near Vr is 
more pronounced at larger W. Hereafter let us set aside the 
unphysical peaky feature near Vr and concentrate on the MFP 
in the subthreshold region (Ag^) and the inversion region 

/ \inv\ 
( A Si?> 

Our key findings in this work with regard to the MFP's 
are: 1) the dimensionless parameter A m /W can be used as a 
good measure of the SRS strength, 2) the MFP's exponentially 
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Fig. 6: (Top) The mean free path A|^ and (Bottom) A™ evaluated at 
Vq - Vr = -0.3 V and 0.4 V versus A m /W, respectively. Both W and 
A m were varied and L m = 1 nm. The data from W = 3, 4, and 5 nm are 
represented by squares, triangles, and diamonds, respectively. The channel 
lengths were fixed to 100 nm for all calculations. The dashed lines represent 
data fit to an exponentially decaying function. 



decrease with A m /W, and 3) Ag^> becomes greater than A^ 
if A m /W > 0.06. 

Firstly, both A|^ and A™^ exponentially decrease with 
A m /W as shown in Fig. __] where Xsr's, were calculated for 
all the possible combinations of W = 3, 4, and 5 nm and A m 
is 0.1 - 0.4 nm. Ag^ and Ag$ were evaluated at V G - V T = 
—0.3 V and 0.4 V, respectively. Notice that W is restricted to 
below 5 nm here. In particular, A^ for W > 5 nm converges 
to the almost the same value, as shown in Fig. __] 

The fact that the MFP can be described by a single param- 
eter A m /W, even with different W's, is quite interesting. For 
the size quantization effects associated with channel width, 
such as the degree of wave-function confinement, the number 
of subbands, and the degree of the mixing between subbands, 
are expected to result in the MFP's dependence on W in 
a complicate way. Nevertheless, our calculations show the 
MFP's exponential dependence on A m /W. This should give a 
useful guide for the estimation of Xsr in Si nanowire devices. 

Secondly, we have found that, if the SRS effects become 
sufficiently strong, Xsr increases with the increase of Vg and, 
as the consequence, A™^ > A|^. This is in sharp contrast to 
the usual case where Xsr should decrease with the increase 
of Vg, because more intense scattering is expected to occur at 
higher gate voltages. In our calculations, the latter is observed 
for the cases of the SRS effects being relatively weak. See 
X S r for W = 5 nm (A m = 0.2 nm) in Fig. __]and that for W 
= 6 and 7 nm (A m = 0.3 nm) in Fig. __] respectively. However, 
as W becomes smaller in Fig. __] the usual behavior is not seen 
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Fig. 7: (a) The mean free path \sr versus Vq for W = 3 nm and L = 100 
nm, as A m is changed from 0.2 to 0.4 nm and L m = 1 nm. (b) The ratio 

T.IW. 




(nm) 



(nm) 



Fig. 8: The cross-sectional electron density profiles for low Vq (left column) 
and high Vq (right column), for (a) W = 7 nm, A m = 0.2 nm and (b) W = 
3 nm, A m = 0.2 nm, and (c) W = 3 nm, A m = 0.4 nm. Electron densities 
are normalized by the total charge density of the cross-section. 



any longer and A™^ becomes greater than X 3 S ^. 

We argue that the ratio A^/A^ has to do with the strength 
of the SRS effects and that A m /W is a good measure for the 
strength. Obviously, the SRS effects should become stronger as 
W becomes smaller for fixed A m (case A) or as A m becomes 
bigger for fixed W^case B). For case A, A m is fixed to 0.3 
nm and W is varied as shown in Fig. [5] and we observe that 
Xsr/Xsr becomes bigger than 1 if A m /W > 0.06. For case 
B, W is fixed to 3 nm and A m is varied as shown in Fig. 
|7](a). For A m = 0.2 nm in the figure, Xsr. already starts to 
ascend in the inversion regime, although quite slowly. As A m 
is increased to 0.3 nm and then to 0.4 nm, \$r ascends more 
and more steeply and A™^ becomes as big as about ten times 
of A|^ for A m = 0.4 nm. (Also compare the curves of W = 
5 nm in Fig. |4] and Fig. [5] where A m is changed from 0.2 nm 
to 0.3 nm.) 

If we combine cases A and B and plot Ag^/Ag^ versus 
A m /W in Fig. (b), we observe that the ratio A^/A^ 
generally increases with the increase of A m JW, although 
there are small fluctuations. A™^/A^ crosses the line of 1 
at A m /W ~ 0.06 in the figure. Thus we conclude that the 
characteristics of the SRS effects in Si nanowires change at 
A m /W ~ 0.06. 

That Xg l ft becomes greater than X s s u ^ for the strong SRS 
effects can be qualitatively explained as follows. In Figs. [H] 
(b) and (c), the cross-sectional charge densities for the cases 
of A m = 0.2 and 0.4 nm are compared at low and high gate 
voltages (W = 3 nm). In the charge density profiles at low 
gate voltage, the Si/SiC>2 interfaces that fluctuate in proportion 
to A,„ can be directly observed. In the charge density profiles 
at high gate voltage, however, the degree of the interface 
fluctuations seems to be almost invariant to the increase of 



the surface roughness. The latter is due to the wave function 
confinement effect which is of particular importance for 
nanowires of W < 5 nm, and consequently, the electrons 
are concentrated to the center of the nanowire as seen in the 
figure ||4), lfl6l . Hence the 'effective' surface roughness that 
felt by the electrons is less at high gate voltage than at low 
gate voltage. This effect intensifies as A m is increased or W 
is decreased, resulting in the relationship between X^/X 6 ^ 
and A m /W shown in Fig. [7] (b). On the other hand, for 
nanowires of W > 5 nm, the confinement effect is weak so 
the electrons are pushed to the nanowire surfaces at high 
gate voltage, resulting in the expected behavior that the MFP 
decreases with the increase of the gate voltage as shown in 
Fig. l](a). 



V. Conclusion 

In this work, we have systematically devised the methodol- 
ogy for the calculation of the mean free path in the nanowire 
field effect transistors with finite channel length. The parasitic 
entrance scattering effects were identified and carefully elim- 
inated from the MFP calculations. Non-perturbative approach 
based on the non-equilibrium Green's function method was 
used to address the SRS effects in the nanowire transistors. 
The methodology developed in this work can be applied to 
any scattering mechanism, although only SRS is considered 
here. 

We have found that the SR-limited MFP can be well 
described by a function which depends on a single parameter 
A m /W. In particular, the MFP in the subthreshold regime 
exponentially decreases with the increase of A m /W. For 



nanowires of W < 5 nm, the MFP in the inversion regime 
also exponentially decreases with A m /W. 

We have also found that the dimensionless parameter 
A m /W can be used as a good measure of the SRS strength. 
If A m /W < 0.06, one may expect that the SRS effects 
should be relatively weak, resulting in the usual behavior 
that the MFP decreases with the increase of the gate electric 
field. If A m /W > 0.06, on the other hand, the SRS effects 
can be regarded to be strong, resulting in the extraordinary 
behavior that the MFP increases with the increase of the gate 
electric field. Our results should give a useful guide for the 
estimation of MFP and the strength of the SRS effects in the 
nanowire devices. 
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